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EXPERIMENTAL EVIDENCE IN SUPPORT OF JOULE HEATING 
ASSOCIATED WITH GEOMAGNETIC ACTIVITY 

SUMMARY 


High resolution accelerometer measurements in the altitude region 
140 to 300 km from a satellite in a near-polar orbit during a period of 
extremely high geomagnetic activity indicate Joule heating to be the 
primary source of energy for atmospheric heating associated with geo- 
magnetic activity. This conclusion is supported by the following 
observational evidence: (1) There is an atmospheric response in the 

auroral zone which is nearly simultaneous with the onset of geomagnetic 
activity, with no s ignificant respons e in the equatorial region until 
several hours later; (2) the maximum heating occurs at geographic loca- 
tions near the maximum current of the auroral electrojet; and (3) there 
is evidence of atmospheric waves originating near the auroral zone at 
altitudes where Joule heating would be expected to occur. An analysis 
of atmospheric respons e time to this heating shows time delays are 
apparently independent of altitude but are strongly dependent , upon 
geomagnetic latitude. 


INTRODUCTION 

Variations in the density of the neutral atmosphere associated with 
geomagnetic activity were first discovered by Jacchia [1960] during an 
investigation to establish an improved relationship between the solar 
radio flux at 10.7 cm and accelerations of satellites. He noted that 
the curve of the satellite accelerations agreed excellently with the plot 
of the solar radio flux except for two transient, short-period secondary 
oscillations. These oscillations were exactly coincident with the dates 
of the only two large geomagnetic disturbances that occurred during the 
lifetime of these particular satellites. Jacchia concluded, from these 
results, that corpuscular radiation capable of producing a large geo- 
magnetic storm has a greater, although more transient, effect on atmos- 
pheric density at 200 km than the fluctuations in solar electromagnetic 
radiation associated with the 10.7 cm solar radio flux. Numerous analyses 
of high altitude density data during the past ten years have shown that 
density increases are almost certain to occur simultaneously with, or 
shortly after, enhancement of geomagnetic activity. 


POSSIBLE ATMOSPHERIC HEATING MECHANISMS 


A number of hypotheses have been proposed to explain the mechanism 
by which this heating occurs, Dessler [1959] advanced the theory that 
such heating results from large-amplitude waves that are formed on the 
edge of the earth’s magnetosphere when it interacts with the solar wind. 
He concluded that, although the electric conductivity of the ambient 
atmosphere at high altitudes is sufficiently high that dissipation of 
such waves due to electron-ion and ion-atom collisions is negligible, 
collisions between ions and neutral atoms at lower levels can cause a 
significant amount of heating, Dessler estimated the altitude at which 
the maximum heating from hydromagnetic waves occurs was between 150 and 
200 km. Akasofu [1965b] subsequently calculated the rise in atmospheric 
temperature due to the hydromagnetic waves which are likely to occur 
naturally in the expected frequency range (uj = 1-25 Hz), and found it was 
not greater than 5 degrees. These results strongly indicate that hydro- 
magnetic waves do not contribute significantly to atmospheric heating 
associated with geomagnetic activity. 

Another possible mechanism suggested to explain heating of the 
ionosphere during magnetic disturbances is corpuscular bombardment by 
energetic charged particles, Mclllwain [1960] suggested that electrons 
in the energy range of 1 to 100 kev are absorbed primarily below 120 km 
altitude. Bremss trahlung from these electrons would presumably be 
absorbed principally at lower altitudes. Cole [1963] calculated an upper 
limit of energy input into the atmosphere by this mechanism to be about 
4 x 10" 6 ergs cm" 3 sec" 1 , and concluded that the energy absorbed by 
neutral particles has little effect on their temperature because of the 
high density of the atmosphere below 120 km. 

Cole [1962] applied the results of earlier work by Cowling and 
Piddington [1954] on Joule heating and motion of uniform ionized gas 
to the ionosphere during geomagnetic disturbances to explain the source 
of energy available during geomagnetic storms by this process. Energetic 
particles precipitate into the auroral zones during periods of enhanced 
geomagnetic activity and provide some direct heat to the neutral thermo- 
sphere by conversion of their kinetic energy. In addition, the increase 
of electron density within the auroral zones results in an increase in 
electric currents within the ionospheric auroral electrojet and in assoc- 
iated heating of the thermosphere at high latitudes by Joule heating. 
Current of density (j) causes Joule heating at the rate Q given by 
j 2 /a 3 , where a 3 is the Cowling conductivity. Thus, since j is propor- 
tional to the magnetic disturbance, the power input is proportional to 
the electric current flowing in the ionosphere. Using the most likely 
estimates of j and a, Cole computed energy sources associated with Joule 
heating. He concluded that they are more than sufficient to explain the 
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heat balance of the disturbed ionosphere. Heating by this process occurs 
in the auroral zones almost simultaneously with the enhancement of geo- 
magnetic activity at altitudes between 100 and 200 km. The heat trans- 
port to lower latitudes could be by heat conduction, heat convection, 
or wave mechanisms. The significant conclusions of this pioneering 
study by Cole [1962] were: 

1. Joule heating in the region 100 to 200 km altitude is common 
during geomagnetic disturbances. 

2. Scale heights and temperatures at altitudes above about 100 km 
increase with geomagnetic disturbances. 

3. The energy of the process causing a geomagnetic disturbance is 
in general not measurable by the energy of the geomagnetic 
disturbance. 

4. A horizontal gradient of pressure at heights above 150 km is 
maintained in the auroral zone during geomagnetic disturbances. 

5. Wind speeds in the polar ionosphere increase during geomagnetic 
d is turbances . 


OBSERVATIONAL EVIDENCE OF ATMOSPHERIC HEATING ASSOCIATED WITH 
GEOMAGNETIC ACTIVITY FROM STUDIES OF SATELLITE DECAY 


Extensive studies of satellite drag-deduced density data by Jacchia 
et al. [1966] revealed several important details of density variations 
associated with geomagnetic activity. Some of these details are listed 
below, 

1. Even small variations in geomagnetic activity are reflected in 
atmospheric density observations. The corresponding temperature varia- 
tions appear related in a near-linear fashion with the planetary geomag- 
netic index, a p , during magnetic storms proper, while during magnetically 
quiet periods the relationship is nearly linear with the planetary index, 
Kp, which is the logarithmic counterpart of ap. According to Jacchia et 
al. [1966], the empirical relationships between these quantities can be 
represented to a fair degree of accuracy by the relationship: 


AT = 1?0 a p + 100°[1 - exp (-0.08 a p ) ] 


( 1 ) 


where AT is the change in exospheric temperature. 
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2. Atmospheric perturbations lag behind geomagnetic perturbations 
by several hours, with an average time lag of 6.7 hours (7.2 hours at 
low latitudes, less than 6 hours at high latitudes). 

3. The temperature during a geomagnetic disturbance is, at least 
occasionally, enhanced at high latitudes. 

These authors recognized that transient density fluctuations such 
as those associated with geomagnetic activity are difficult to study by 
the satellite drag-analysis method, because the time and spatial resolu- 
tions (about 5 hours or 3 orbits) are so large. Since the quantity to 
be determined, the acceleration of the satellite* s mean motion, is the 
second derivative of the mean anomaly, this task becomes increasingly 
difficult. 

Another very important limitation in the technique of deducing 
density values from satellite drag data is the fact that the change in 
the orbital period is a measure of the integrated effect of air drag in 
the neighborhood of perigee. This angular resolution is a function of 
the orbit eccentricity. For example, King-Hele [1966] has calculated 
that for an orbit with an eccentricity of .2 and a perigee altitude of 
200-250 km, the drag is appreciable over an arc of 20° on each side of 
perigee. Since this drag occurs at heights up to one scale height, H 
(about 25 km at altitudes of 200-250 km), above perigee, the resulting 
density value represents, assuming a time resolution of .2 day, an 
average through a region up to 25 km in height over a lateral geographic 
distance up to 3400 km. Thus, it is not surprising that analyses of 
these data lead to different conclusions concerning time delays of 
density increases following geomagnetic storms. 

To identify the energy source and the mechanism by which the atmos- 
phere is heated following geomagnetic storms, the intensity of the heat- 
ing and the reaction time of the atmosphere are two parameters that are 
especially needed. From an analysis of three high inclination satellites, 
Injun 3, Explorer 19, and Explorer 24, Roemer [1969] made some very 
significant contributions concerning the intensity of heating of the 
atmosphere for a given level of geomagnetic disturbance. Although con- 
siderable scatter exists in his data, his results indicate that the 
heating level (AT/z^a p ) at auroral latitudes is about twice as large as 
that at the equator. In a study of the reaction time of the atmosphere 
to heating associated with geomagnetic activity, Roemer compared times 
at which the rate of change of satellite periods reached their maximum 
with the times of occurrences of the peak of geomagnetic variations for 
about 100 geomagnetic events. This analysis indicated a mean time delay 
time. At = 5.3 ± 0.4 (s.d.) hours. These average results concerning time 
delays are in fair agreement with those obtained by Jacchia et al . [1966], 

who analyzed 80 events from four high inclination satellites and found an 
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average lag time in the low to medium latitude region (0 < 55°) of 7.2 
hours, and an average time lag in the high latitude region (0 > 55°) of 
5.8 hours. There was considerable scatter in these results, with time 
delays ranging from zero to more than 14 hours. 

DeVries et al. [1967] analyzed high-resolution density data deduced 
from precise radar tracking data from eleven low-altitude, attitude- 
stabilized Agena satellites and discovered delay times of density 
increases associated with geomagnetic activity from near-instantaneous 
in the auroral regions (0 = 75 degrees) to about 17 hours at lower lati- 
tudes (0 = 30 degrees). In this study, 27 events during which the geo- 
magnetic planetary index, a p , ranged from 12 to 154 units were analyzed. 
One recognized deficiency in these data was the relatively small eccen- 
tricities (.013 - .027) of the orbits of the satellites under considera- 
tion. This deficiency increases the angle about perigee through which 
the drag occurs. Jacchia et al. [1966] noted, and DeVries agrees, that 
the low eccentricities make it difficult to assign a latitude to a drag 
observation density. 

Incorrect conclusions concerning time delays between geomagnetic 
disturbances and associated atmospheric heating can result from improper 
smoothing of geomagnetic data. Jacobs [1967] analyzed the response time 
of the atmosphere to geomagnetic disturbances in the 165-200 km region 
of the atmosphere by assuming that the geomagnetic indices could be 
represented by 12-hour running means of the three-hourly geomagnetic 
index, K p , appropriate for a time delay for the atmospheric response of 
.23 day. Lew [1969] later reanalyzed the same data and found that the 
actual response time Jacobs had assumed was .48 day, because of the 
manner in which he chose to assign the K p averages. By averaging the 
K p values over a period of 12 hours (4 data points) and assigning the 
average value to the last time point, Jacobs had inadvertently intro- 
duced a time displacement error of about six hours. _In s itu measure- 
ments of neutral atmospheric composition and density, which offer time 
resolutions on the order of seconds as compared to hours for drag 
derived data, have recently become available. Taeusch et al. [1970] 
analyzed neutral atmospheric composition variations obtained from a 
quadrupole mass analyzer during a period in which several geomagnetic 
disturbances occurred. They found the response of the atmosphere to 
this activity to be localized in the auroral regions of the earth, with 
a response time of less than one hour. The short response time at the 
high latitudes is in good agreement with that found by DeVries [1967] 
at comparable latitudes. 
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LOGACS - AN ORBITAL ACCELEROMETER CALIBRATION EXPERIMENT 


One logical means of measuring satellite accelerations caused by 
atmospheric variations is by an on-board accelerometer. A major diffi- 
culty in this measurement method is the difficulty in calibrating 
accurately the accelerometer for the low accelerations in an orbital 
environment. To overcome this difficulty, in 1967 the Air Force formu- 
lated an experiment to calibrate an extremely sensitive accelerometer in 
the orbital environment. This experiment is known as LOGACS (LOw-G 
Accelerometer Calibration System). 

The accelerometer selected for the LOGACS flight was the Bell MESA 
(Miniature Electrostatic Accelerometer). The MESA is a single-axis, 
electrostatically pulse-rebalanced accelerometer designed so that both 
the suspension and rebalance forces can be changed, or scaled, by chang- 
ing the voltage levels. For this experiment, the Bell MESA was mounted 
in the aft region of a standard Agena vehicle. The accelerometer was 
placed on a turntable device so that it could be rotated at two pre- 
cisely known angular rates or held in two known fixed positions. The 
Agena vehicle was attitude controlled such that the plane of the turn- 
table was always held in the plane of the local horizon. In its fixed 
positions, the sensitive axis of the MESA was aligned in the orbital 
plane (one position in the direction of motion and the other opposite 
to the direction of motion). Comparison of these measurements allowed 
accurate determination of the instrument bias. The instrument scale 
factor was determined by rotating the accelerometer at two distinct 
rotation rates, approximately 1/2 and 1 rpm. A comparison of measure- 
ments obtained in each of these modes permitted accurate determination 
of the scale factor. Fotou [1968] has prepared a detailed description 
of the instrumentation, data processing methods, results, and other 
pertinent details concerning the experiment. 

The LOGACS experiment was placed into a near-polar orbit on May 22, 
1967. During the 18th revolution one day later, the orbit was adjusted 
such that its apogee altitude was raised to increase the orbital lifetime. 
The orbital parameters at injection and after the orbit adjustment are 
given in Table 1. 
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Table 1. Orbital Parameters for LOGACS after Orbit Adjust 



Injection 

Rev. 18.4 

Date 

May 22, 1967 

May 23, 1967 

Time 

1839 GMT 

2124 GMT 

Height of perigee 

148 km 

147 km 

Height of apogee 

357 km 

403 km 

Inclination 

91.5 deg 

91.5 deg 

Eccentricity 

0.015 

0.018 

Period 

89.4 min 

89.8 min 

Latitude of perigee 

43.3 deg N 

40.7 deg N 

TCpCTTT TO r\T? 

l\b JUljllJ W i- 

ANALYSES OF LOGACS 

DATA 


The LOGACS experiment, which was a complete success, provided 
approximately 100 hours of acceleration data from which the null bias, 
accelerometer scale factor, and atmospheric drag on the satellite were 
calculated. Bruce [1968] estimates the actual errors in the accelera- 
tions measured during the LOGACS experiment are less than one percent 
and the density values deduced from the acceleration measurements are 
accurate to within plus or minus 10 percent. The major source of 
error is the uncertainty associated with the drag coefficient. 


Unexpected geophysical events during the LOGACS experiment were the 
occurrences of two class 2B solar flares and one class 3 b solar flare 
on May 23, 1967, followed by greatly enhanced geomagnetic activity. 
Geomagnetic activity was fairly quiet the first portion of the flight, 
but a large geomagnetic storm, with a geomagnetic index, a p , value of 
256 units, during revolution 44 was followed about nine hours later 
by a great geomagnetic storm, with geomagnetic a p value of 400 for two 
successive 3-hour reporting periods. The values of the 10.7 cm solar 
radio flux and the geomagnetic planetary index, a p , during the LOGACS 
flight are shown in Figure 1, along with the LOGACS orbit numbers 
during which useful density data were derived. 

1. Comparison of Drag-Deduced and Accelerometer Density Data . The 
LOGACS vehicle, which also carried an S-band transponder, was tracked 
accurately by the pulse radar network of the U. S. Air Force Satellite 
Control Facility. Drag-deduced data during the experiment were computed 
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for comparison with the accelerometer data. The drag-deduced densities 
and the accelerometer-deduced densities are shown in Figure 2. This 
graph reveals some interesting differences in the density values obtained 
by the two methods, such as a difference of about 3.6 hours between the 
peak density values detected by the accelerometer and those computed 
from the drag-deduced data during the two periods of enhanced geomagnetic 
activity and differences of density values up to 15 percent between the 
two methods of measurement. 

These apparent discrepancies are not surprising, inasmuch as the 
drag-deduced data are smoothed over periods of several hours and over 
large distances from the perigee position, while the accelerometer data 
are near-instantaneous measurements. As would be expected, the accelerom- 
eter measurements are much more sensitive to atmosphere variations than 
drag-deduced measurements because of their greatly improved time and 
spatial resolution. 

2 . Latitudinal Response of the Atmosphere to Corpuscular Heating . 

The orbit geometry of the LOGACS experiment provided a unique oppor- 
tunity to study the atmospheric response associated with geomagnetic 
activity as a function of geomagnetic latitude. The altitudes of the 
LOGACS satellite for a typical orbit are shown in Figure 3. During a 
typical orbit, the satellite was at an altitude of about 300 km at the 
ascending node, descending to 200 km over the north pole and to about 
145 km at perigee, then ascending to about 170 km over the equator 
and to about 300 km over the south pole. 

The highest resolution density data previously available from satel- 
lite orbit decay was on the order of about three hours, or two orbits; 
LOGACS provided data with a time resolution of 20 seconds -- an improve- 
ment in time resolution of a factor of more than 500. 

During the first 40 orbits of the LOGACS experiment, the geomagnetic 
activity was low, with the values of the geomagnetic planetary index, ap, 
not exceeding 32 units. Fortunately, they were seven units or less for 
nine hours preceding the two large geomagnetic storms that occurred 
during the flight. The three-hourly a p values increased from four units 
at orbit 41 to 56 units three hours later and to 236 units six hours 
later. A decrease of a p values to 154 units followed, with an increase 
to 400 units during orbit 51. 

For this analysis of density increases associated with the enhanced 
geomagnetic activity during the LOGACS experiment, density data subsequent 
to orbit 41 were normalized to the geomagnetically quiet conditions occur- 
ring during orbit 41. The ratios of the normalized densities as a func- 
tion of geomagnetic latitude were then plotted for all orbits for which 
data were available. 
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The satellite positions corresponding to the left sides of Figures 
4-9 were near the ascending node of the orbit, or the points at which 
the satellite crossed the equatorial plane in a northerly direction. 

At these positions, the satellite would be at an altitude of about 
300 km. The first crossing of the northern hemisphere auroral zone 
would occur on the night sector, near 63-70 degrees geomagnetic lati- 
tude, at an altitude of about 190 km. After crossing the north pole, 
the satellite would cross the auroral zone on the sunlit side of the 
hemisphere. The satellite was in a near sun-synchronous orbit, with 
a local time at perigee of about 10:30. Ratios of the density at each 
revolution to that occurring during revolution 41 are shown in the 
upper part of the ordinates of Figures 4-9. The altitudes of the 
satellite are shown in the bottom part of the ordinates, with the geo- 
magnetic latitudes shown in the abscissas of Figures 4-9. 

Figure 4 shows the graph of the normalized density versus geo- 
magnetic latitude during revolution 43, at which time the value of the 
a p index had increased to 56 units. In general, the density values had 
not changed more than a few percent at geomagnetic latitudes less than 
50 degrees. However, the density at a geomagnetic latitude of about 
75 degrees (on the dark side of the northern hemisphere) had increased 
about 30 percent, while that at the sunlit side of the northern hemi- 
sphere at the same latitude had decreased about 20 percent. This is 
a rather surprising result, since the most widely accepted atmospheric 
models such as the Jacchia [1970] and CIRA [1965] would indicate global 
heating about six hours subsequent to the enhanced geomagnetic activity. 

Figure 5 shows the neutral density enhancement during revolution 
46, approximately 1 1/2 hours subsequent to the peak of the first large 
geomagnetic storm, during which the a p value reached 256 units. Again, 
the highest relative increase of density during this orbit occurred at 
about 70 degrees geomagnetic latitude on the night side of the northern 
hemisphere. A very small (about 5 percent) increase of density occurred 
at 70 degrees geomagnetic latitude on the sunlit side; however, signifi- 
cant increases occurred at higher geomagnetic latitudes in the southern 
hemisphere. This indicates that the immediate source of heating is near 
the auroral zone on the night side of the northern hemisphere and near 
the auroral zone in the southern hemisphere. 

Figure 6 shows the normalized density values during revolution 47, 
about three hours subsequent to the peak of the first geomagnetic storm. 

In general, density values at lower latitudes were starting to increase 
(by about 20 percent), while those in the auroral regions were decreasing. 

Figure 7 shows the normalized density ratios at about 41/2 hours 
subsequent to the first geomagnetic storm. By this time, density values 
near the equator had increased about 50 percent, while those near the 
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auroral regions on the night side of the northern hemisphere had decreased 
to values near those at geomagnetical ly quiet times. One very interest- 
ing feature illustrated in Figure 7 is the decrease of density over the 
polar cap and in the auroral regions to values 20 to 30 percent below 
those occurring before the enhanced geomagnetic activity. 

About nine hours after the peak geomagnetic index values for the 
first geomagnetic storm had been reached, one of the largest disturbances 
of solar cycle 19 occurred, with values of the geomagnetic index, ap, 
reaching 400 units for two consecutive periods. Figure 8 shows ratios 
of densities at the beginning of this high geomagnetic activity to those 
occurring during the geomagnetically quiet period 18 hours previously. 
Again, the significant density increases were localized near the geo- 
magnetic latitude of 65 degrees on the dark side of the northern hemi- 
sphere and near 60 degrees in the observations available from the 
southern hemisphere. In both cases, values of density increased by 
more than a factor of two, almost simultaneously with the increase of 
geomagnetic activity. An examination of the satellite altitudes at which 
the highest density ratios were evident reveals no noticeable altitude 
dependence, since the satellite was at an altitude of about 190 km during 
the northern hemisphere traverse of the night auroral zone and at an 
altitude of about 265 km during the period of enhanced density values 
in the southern hemisphere auroral zone. Again, there was a decrease of 
density to below geomagnetically quiet values at a geomagnetic latitude 
of about 65 degrees on the sunlit side of the northern hemisphere. 

Density values near low altitudes, at this time, had not changed signifi- 
cantly from slightly enhanced values associated with the previous geo- 
magnetic activity. 

The latitudinal variations of density changed considerably within 
the next 4 1/2 hours, as shown in Figure 9. Although the a p index during 
this orbit had remained at 400 units, the atmospheric density decreased 
to near-normal values in the auroral region, but had increased by more 
than a factor of about two in equatorial regions. Peak density ratios 
had thus moved from 65 degrees latitude to 20-30 degrees latitude during 
the 4 1/2 hour period, although there had been no decrease in geomagnetic 
activity. At the auroral latitudes on the sunlit side of the northern 
hemisphere, the density had decreased to about 75 percent of that occurring 
before the onset of geomagnetic activity. 

3 . Dependence on Time Delays for Atmospheric Heating on Geomagnetic 
Latitude and Altitude . During the 36-hour period for which density were 
available from the LOGACS experiment, following the onset of enhanced 
geomagnetic activity, density values were obtainable from all or portions 
of 20 orbits. These data presented a unique opportunity to evaluate 
effects of geomagnetic latitude and altitude on the time required for the 
neutral atmosphere to respond to the heating associated with geomagnetic 
activity. Density data normalized to the geomagnetically quiet period 
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preceding the high geomagnetic activity, as described in the previous 
section, were tabulated at intervals of 10 degrees geomagnetic latitude 
for each available orbit* Time cross sections for each 10 degrees of 
geomagnetic latitude were then prepared for the entire period. Since 
the orbital periods were about 90 minutes, this was the minimum time 
resolution attainable. The longest time resolution during the period 
under consideration was three orbits, or about 270 minutes. 

Time cross sections of normalized density values at intervals of 
10 degrees geomagnetic latitude are shown in Figures 11 and 12. Corre- 
sponding values of the geomagnetic index, ap, are shown at the bottom 
of the figures. Time delays between the peak of the a p values for each 
of the two large geomagnetic events that occurred during this period 
(assumed at the middle of the 3-hourly a p reporting periods) and the 
peak values of the normalized density associated with each event are 
shown for all cases in which adequate data were available to make this 
evaluation. 

An analysis of these time delays strongly supports the earlier 
conclusions by DeVries et al. [1967] that time delays between the onset 
of geomagnetic activity and the subsequent enhancement of atmospheric 
heating are latitude-dependent. For example, at geomagnetic latitudes 
greater than 60 degrees, these time delays were, in every case, not 
more than 1 1/2 hours; at geomagnetic latitudes less than 30 degrees, 
they were, in every case, not less than 4 1/2 hours. Between these 
latitudes, the time delays were about in proportion to the decrease in 
geomagnetic latitudes . 

The time differences between the peak of the geomagnetic disturb- 
ances and the associated atmospheric heating for the 30 cases for which 
data were available are shown in Figure 12. Corresponding altitudes are 
shown in parenthesis. It can be seen from Figure 12 and Figures 5 and 
8 that the time delays are about independent of altitude but are strongly 
dependent upon geomagnetic latitude. 

The serial correlation coefficient of geomagnetic latitudes with 
lag times for the 30 cases is .834, which is less than the .001 level 
of statistical significance. In other words, the probability of 
exceeding this value by chance is less than one in one thousand. The 
line of least squares fit, as shown in Figure 21, is based upon the 
relationship: 


70.4 - Geomagnetic Latitude 
9.4 


( 2 ) 


where T = lag time in hours. 
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4. The Similarity Between the Locations of the Auroral Electrojet 
and Spatial Density Variations Detected During the LQGACS Experiment . 
Auroral electrojets are intense electric currents in the auroral zones 
associated with geomagnetic disturbances. Clumen et al. [1969] summarize 
the general features of the auroral electrojet to consist of a height of 
the jet axis varying from 100 to 150 km, a length from 2000 to 5000 km, 
north-south thickness from 500 to 680 km, vertical half-width of 30 km, 
and a duration of 30-60 minutes. Akasofu et al. [1965a] discuss the 
morphology, geometry, and current system of the auroral electrojet. If 
the mechanism for atmospheric heating associated with geomagnetic dis- 
turbances is Joule heating, according to Cole [1962] the magnitude of 
such heating will be a function of the electrical current in the iono- 
sphere with major concentration in the auroral latitudes. Accepting 
this hypothesis and a model current system for an intense polar magnetic 
substorm similar to one depicted by Akasofu et al. [1965a], the most 
intense heating should occur in the auroral zone of the midnight sector 
of the hemisphere. 

A proposed model current system for an intense magnetic substorm by 
Akasofu et al. [1965a], and the corresponding traverse of the LOGACS 
experiment through the system are shown in Figure 13. LOGACS experiment 
revolutions 46 and 51 occurred shortly after the onsets of two intense 
geomagnetic disturbances. Figures 5 and 8 show the immediate atmospheric 
density enhancements associated with this activity. The portion of the 
orbit in Figures 5 and 8 which corresponds to the track shown in Figure 13 
includes the region from geomagnetic latitudes of plus 50 degrees to plus 
90 degrees, then back to plus 50 degrees. Figures 5 and 8 with the 
orbital track shown in Figure 13 show an interesting pattern of density 
enhancements which resembles to a remarkable extent that which would be 
expected with the current system in Akasofu ? s model. In both cases, the 
most intense enhancement of density occurs about coincident with the 
most intense portion of the electrojet current system penetrated by the 
satellite (Point A, at about 70 degrees geomagnetic latitude at an alti- 
tude of 190 km in Figure 5 and at about 75 degrees geomagnetic latitude 
and an altitude of about 190 km in Figure 8). The portions of the orbits 
shown in Figures 5 and 8 from point B to C (75 to 50 degrees geomagnetic 
latitude on the sunlit side of the hemisphere) show no significant enhance- 
ment of density even though they are in close proximity to the auroral 
oval. An examination of this region in Akasofu 1 s model (Figure 13) shows 
there is no enhancement of the current system in this region; therefore, 
Joule heating would not provide the mechanism for heating on this portion 
of the orbits. 

Unexpected anomalies, which are shown between points D and E in Fig- 
ure 7 (about 4 1/2 hours subsequent to the commencement of the first geo- 
magnetic storm) and Figure 9 (about 41/2 hours subsequent to the 
commencement of the second geomagnetic storm), occur in regions in which 
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the density decreased to values as low as 70 percent of pre-storm values 
The geographic regions in which this phenomenon occurs appear to be 
within the auroral oval. The apparent density depletion may result from 
the circulation and wind patterns associated with the high geomagnetic 
activity occurring a few hours previously. 

5. Evidence of Wave Propagation from LOGACS Data . In a study to 
evaluate the capability of various atmospheric models to depict atmos- 
pheric density for geomagnetic conditions during the LOGACS flight, 
cross sections of the ratios of the actual density to the density 
depicted by various models as a function of geomagnetic latitude were 
prepared. Figure 14 shows such a cross section during the period sub- 
sequent to the second geomagnetic storm. These data indicate evidence 
of atmospheric waves. General features of these waves include wave- 
lengths on the order of 450 km and amplitudes of the ratio of the actual 
density to that depicted by the model increasing with altitude. It 
seems plausible that the source of this apparent wave motion is Joule 
heating associated with the geomagnetic storm. Such a wave would most 
likely propagate radially from the auroral electrojet. The satellite 
descended in altitude as it moved north through the wave region and the 
wavelength is a combination of vertical and horizontal variations. It 
is noteworthy that the amplitude of the wave decreased as the altitude 
decreased. As the wave propagates from the source it will tend to con- 
serve kinetic energy and the amplitude will decrease on the portion of 
the wave front that is propagating into the more dense layers and 
increase on the portion of the wave that is propagating upward into the 
less dense layers. Theoretical studies on the phenomena are being con- 
ducted in collaboration with Dr, George H. Fichtl of the NASA/MSFC Aero- 
space Environment Division. 

6 . Inadequacy of Currently Available Atmosphere Models to Depict 
Density Variations Associated with Large Geomagnetic Disturbances . The 
most widely accepted model of the atmosphere, the Jacchia Model [1970], 
approximates temperature increases (with a lag of 6 hours) associated 
with geomagnetic activity by equation (1). Actual density increases are 
then computed from a model. The maximum temperature correction by this 
method for the largest geomagnetic storm that can be measured is about 
500 degrees. While this temperature correction appears to be adequate 
to increase the densities sufficiently at altitudes above about 250 km, 
it is impossible to accurately model density increases associated with 
large geomagnetic disturbances at altitudes below about 200 km by this 
technique. Figure 15 shows the extremes of density measured by the 
LOGACS experiment between altitudes of 145 and 290 km, compared to those 
depicted by the most recent Jacchia Model [1970], assuming exospheric 
temperatures of 1000, 1500, and 2000 degrees. Figure 15 shows that maxi 
mum density values observed during the LOGACS experiment when the ap 
value was 400 units are higher than those depicted by the Jacchia 
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2000-degree exospheric temperature model at all altitudes up to 290 km, 
with the percentage differences being larger at lower altitudes. In 
other words, all observed density values when a p was 400 units lie out- 
side the limits of those depicted by the Jacchia Model. 

The method used by Jacchia and others to study the magnitude of 
atmospheric variations is to convert the density variations, as depicted 
by satellite drag, to temperature variations in order to relate events 
at different altitudes to each other. This transformation is only as 
accurate as the atmospheric model being used for the transformation. 

One of the major shortcomings of the currently available atmospheric 
models, including Jacchia 1 s [1970], is that they assume static equili- 
brium conditions which we know do not exist during the short-period 
atmospheric perturbations occurring during a geomagnetic storm. 

The second major problem in the modeling of atmospheric variations 
associated with geomagnetic activity is depicting the reaction time of 
the atmosphere. Current models are designed on the premise that the 
atmosphere reacts on a global basis with a time lag of 6 hours subsequent 
to the geomagnetic activity. Analyses of LOGACS data strongly indicate 
that neither assumption is entirely correct; i.e., the atmosphere responds 
nearly simultaneously with the onset of geomagnetic activity near the 
auroral zones, with this energy being propagated to lower latitudes by 
wave mechanisms, and/or conduction and convection with a time delay of 
about six hours near the equator. 

Based upon these data, the modeling of atmospheric variations assoc- 
iated with geomagnetic disturbances can be improved by incorporating a 
variable time delay (near-simultaneous at the auroral region and increas- 
ing to six hours at the equator) and by directly modeling the density 
perturbations as a function of geomagnetic activity near the altitude 
region where the heating occurs (below 140 km). Research on this problem 
is underway at the MSFC Aerospace Environment Division. Additional data 
from other accelerometer flights will be quite valuable for this effort. 


CONCLUSIONS 


The LOGACS experiment probably furnished more evidence from which to 
evaluate the mechanism responsible for atmospheric heating associated with 
geomagnetic activity than any other data available up to the present time. 
Analyses of these data strongly indicate Joule heating as the principal 
mechanism for such atmospheric response. This conclusion is supported by 
the following observational evidence from the LOGACS experiment. 

1. There is a response of the atmosphere in the auroral regions 
nearly coincident with the onset of geomagnetic activity, with no 
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significant response evident at lower latitudes until several hours 
later . 


2. There is a high correlation between the times of peak geo- 
magnetic activity and the geomagnetic latitudes at which the maximum 
heating of the neutral atmospheric heating occurred. This relation- 
ship indicates that heating associated with the geomagnetic activity 
originated in the auroral region and propagated systematically to lower 
latitudes . 

3. The maximum atmospheric heating occurred at geographic loca- 
tions near the expected maximum current of the auroral electrojet. 

This would be expected with Joule heating, since Q, or energy from 
such heating, is directly proportional to the square of the auroral 
electrojet current. 

4. Strong evidence of atmospheric wave propagation was present 
only during short periods subsequent to enhanced geomagnetic activity. 
These waves appear to originate in the auroral region at altitudes 
commensurate with a source associated with Joule heating. 

5. A decrease in atmospheric density in the auroral regions 
occurred about 4 to 5 hours after commencement of the high geomagnetic 
activity, during a period in which the geomagnetic activity remained 
highly enhanced. This indicates atmospheric circulation associated 
with extremely high pressure gradients and winds which would be 
expected with Joule heating. 

Presently available atmospheric models do not depict adequately 
the atmospheric response to geomagnetic activity, particularly at low 
(below about 200 km) altitudes. This is not surprising, since these 
models are based upon the assumption of equilibrium conditions, which 
are known not to exist during transient atmospheric disturbances. 
Atmospheric models can be improved by incorporating an atmospheric 
response near the altitude at which the heating occurs, with a latitude- 
dependent magnitude and a time delay. 
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Figure 2. Comparison of density deduced from tracking data during the 
LOGACS flight with density deduced from accelerometer data. 
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Figure 3. LOGACS Orbit Geometry 
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Figure 5. Normalized density values versus geomagnetic latitudes during revolution 46, 
about 1 1/2 hours after a n = 256. 
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Figure 6. Normalized density values versus geomagnetic latitudes during revolution 47, about 
3 hours after a n - 256. 
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Figure 7. Normalized density values versus geomagnetic latitudes during revolution 48, about 
41/2 hours after a D = 256. 
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Figure 10. Time cross-sections of normalized density values at 10 degrees 
(0=0- ±30 degrees geomagnetic latitude). 
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Figure 11 # Time cross sections of normalized density values at 10 degrees of latitude intervals 
(O = ± 40 - ± 80 degrees geomagnetic latitude). 











Figure 13. Proposed model current system for intense polar magnetic 
storms, with the corresponding LOGACS track (after 
Akasofu [1965a]). 
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Figure 14. Evidence from LOGACS data of atmospheric wave propagation. 
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Figure 15. Extreme values of density depicted by the Jacchia 70 atmos- 
phere compared to those observed by the LOGACS experiment. 
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